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We present an initial analysis of AKARI 18 /im-selected galaxies using all 9 photometric bands at 2 
- 24yLtm available in the InfraRed Camera (IRC), in order to demonstrate new capabilities of AKARI 
' cosmological surveys with this unprecedented wavelength coverage at mid-infrared bands. We detected 72 

, sources at 18 //m in an area of 50.2 arcmin^ in the AKARI/IRC monitor field towards the North Ecliptic 

Pole (NEP). From this sample, 25 galaxies with probable redshifts z ^ 0.5 are selected with a single colour 
cut (iV2-iV3 > 0.1) for a detailed SED analysis with ground-based BVRi'z'JK data. Using an SED 
radiative transfer model of starbursts covering the wavelength range UV ~- submm, we derive photometric 
5— i ' redshifts from the optical-MIR SEDs of 18 /^m-selected galaxies. From the best-fit SED models, wc show 

that the IRC all-band photometry is capable of tracing the steep rise in flux at the blue side of the PAH 
6.2 fim emission feature. This indicates that the IRC all-band photometry is useful to constrain the rcdshift 
of infrared galaxies, specifically for dusty galaxies with a less prominent 4000 A break. Also, we find that 
the flux dip between the PAH 7.7 and 11.2 /im emission feature is recognizable in the observed SEDs of 
galaxies at z ~ 1. By using such a colour anomaly due to the PAH and silicate absorption features, unique 
samples of ULIRGs at z ~ 1, 'silicate-break' galaxies, can be constructed from large cosmological surveys of 
AKARI towards the NEP, i.e. the NEP-Dcep and NEP- Wide survey. For AGN candidates selected with an 
IRC colour-colour diagram (iV2 — A^3 vs. N3 — S7), we find excess emission in more than two MIR bands, 
compared to the best-fit starburst SED model, which suggests the presence of hot dust around the central 
black hole. This pilot study suggests the possibility of detecting many interesting galaxy properties in the 
NEP-Deep and Wide surveys, such as a systematic difference in SEDs between high- and low-z ULIRGs, 
and a large variation of the PAH inter-band strength ratio in galaxies at high redshifts. 

Key words: infrared: galaxies — galaxies: starburst — galaxies: active — galaxies : evolution 



2 



T. Takagi et al. 



1. Introduction 

The universe is a much more luminous place than 
seen in the optical, so it has been revealed since the 
launch of the first infrared astronomical satellite, IRAS 
(Neugebauer et al., 1984). Subsequent infrared space mis- 
sions, i.e. ISO (Kessler et al., 1996) and Spitzer (Werner 
et al., 2004) have led to enormous advances in our knowl- 
edge on star-formation activity, specifically in the distant 
universe (e.g. Genzel & Cesarsky, 2000; Franceschini et al., 
2006, and references therein). The last mission, AKARI 
was launched on 22nd February 2006 JST (Murakami 
et al., 2007) and is functioning very well in orbit at the 
time of writing. 

The deepest extragalactic surveys for infrared luminous 
galaxies have been carried out at 15 ^m with ISO and 
24 /im with Spitzer, due to the fact that the confusion of 
point sources and of background cirrus emission is more 
serious in longer wavelengths. A key result from ISO sur- 
veys may be summarised in the differential number counts 
at 15/im (see e.g. Figure 14 in Genzel & Cesarsky, 2000), 
which have a characteristic bump at around 300 i-iJy. With 
these number counts, pure density evolution models are 
clearly rejected, and strong luminosity evolution with 
~ (1 -|- z)'^ is suggested at < z < 1 (Xu, 2000; Pearson, 
2001; Franceschini et al., 2001; Chary & Elbaz, 2001). 
Recently, this strong luminosity evolution at < z < 1 
has been confirmed by Spitzer with much greater statis- 
tical significance (Perez-Gonzalez et al., 2005; Babbedge 
et al., 2006; Caputi et al., 2007). 

Multi- wavelength surveys at both 15 and 24 /im are not 
just redundant observations, but reveal the problems in 
our understanding of distant infrared galaxies. Although 
the differential number counts at 24 fiia are found to be 
similar in shape to the 15 fim counts, it is not straightfor- 
ward to explain the 24 fim counts with the phenomenolog- 
ical models of galaxy counts constrained to fit the multi- 
wavelength counts prior to the Spitzer results. The peak 
in the 24 /im differential counts is found at 200 - 300 /t Jy 
(in dN/dS x 5^'^) which is fainter than expected from 
models of galaxy counts. The resolution of this problem 
has not yet been achieved (see however Lagache et al., 
2004; Chary et al., 2004; Pearson, 2005). Since at 15 and 
24 /im we observe the PAH emission features of galaxies 
at z =1 - 3, we need SED templates with realistic PAH 
emission features as a function of luminosity for galaxies 
at z > 1. Unfortunately, such SED templates are not yet 
available. 

With the InfraRed Spectrograph (IRS) onboard Spitzer, 
MIR features of high-z galaxies, such as the PAH emission 
feature and silicate absorption feature, become accessible 
(Houck et al., 2005; Lutz et al., 2005; Yan et al, 2005; 
Huang et al., 2006; Weedman et al., 2006; Menendez- 
Delmestre et al., 2007; Teplitz et al., 2007), although the 
observed samples are still very limited. The highest red- 
shift PAH detection has been reported in an infrared- 
luminous Lyman break galaxy at z = 3.01 (Huang et al., 
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2006) , which is likely to be a hypcrluminous infrared 
galaxy. Currently, the prime high-z targets for IRS 
seem to be submm galaxies (Lutz et al., 2005; Mcncndez- 
Delmestre et al., 2007; Pope et al., 2006). It is rather sur- 
prising that PAH emission features are detected in most 
submm galaxies, given that they are typically hyper lu- 
minous (L/_R > 10^'^Lq) galaxies which are in the local 
universe generally AGN-dominated galaxies. There is al- 
ready a hint that the MIR properties of high-z galaxies 
are systematically different to those in the local universe. 
More systematic studies on this issue are clearly needed. 

With Spitzer surveys, the majority of Hyper- and Ultra- 
Luminous InfraRed Galaxy (Hy/ULIRGs)^ targets at 
high redshifts may be derived from the 24 /im source cata- 
logues. This may result in incomplete samples of ULIRGs 
specifically at z 1.5, at which ULIRGs with strong sil- 
icate absorption become very faint in the 24 /im band 
(Takagi & Pearson, 2005; Armus et al., 2007). Caputi 
et al. (2006) obtained a bimodal redshift distribution of 
24 /im sources as a result of the deficit of galaxies at z ~ 1.5 
(see also Franceschini et al., 2006). More comprehensive 
wavelength coverage at MIR bands should be useful to 
produce complete samples of ULIRGs at 1 < z < 2. We 
point out below that the AKARI multi-wavelength sur- 
veys are best suited for this purpose, thanks to their su- 
perb wavelength coverage in the NIR ~ MIR. 

The most unique scientific programme of AKARI is the 
all-sky survey, the first major update of IRAS. Besides 
this, another uniqueness of AKARI can be found in sur- 
veys of selected areas with unprecedented wavelength cov- 
erage in the NIR ~ MIR, with 9 photometric bands from 2 
to 24 /im using the InfraRed Camera (IRC; Onaka et al., 

2007) . With the IRC, two main extragalactic surveys have 
been designed and conducted around the North Ecliptic 
Pole (NEP), a deep ~ 0.5 deg^ survey (the NEP-Deep sur- 
vey) and a shallow 6.2 deg^ survey (the NEP- Wide sur- 
vey) . See Matsuhara et al. (2006) for details on the survey 
design of the NEP-Deep and NEP- Wide surveys. 

With this comprehensive wavelength coverage, one clear 
benefit for extragalactic MIR surveys is that a more re- 
liable selection for high-z ULIRGs becomes possible, as 
noted above. In a flux limited sample with Spitzer MIPS 
24 /im, it is likely that some fraction of ULIRGs at z ~ 1.5 
are missed, owing to redshifted silicate absorption fea- 
tures. With IRC, such ULIRGs would be detected at 
15 and/or 18 //m. These ULIRGs may be called 24 //m 
dropouts or 'silicate-break' galaxies (Takagi & Pearson, 
2005). This capability would be very useful to make a 
complete sample of MIR-bright ULIRGs for follow-up ob- 
servations with MIR spectroscopy. Furthermore, the re- 
sulting well-sampled SEDs of MIR-selected galaxies allow 
us to study the nature of distant infrared galaxies in detail 
with the help of SED models, including the photometric 
estimation of redshifts from MIR features. 

As a pilot study for a full analysis of the NEP survey, 
here we report the first insight on AKARI 18 /im-selected 

^ [Hy,U]LIRGs defined with logLm/LQ =[> 13, 12 - 13], where 
LiR is the luminosity in the wavelength rage from 8 to 1000 fim. 
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distant infrared galaxies seen by all 9 photometric bands 
available in the IRC. Lee et al. (2007) provide similar 
study for an AKARI 11 /im-sclected sample with 6 pho- 
tometric bands between 2 and II fim. Matsuhara et al. 
(2007) report the optical properties of an AKARI 15 /im- 
selected sample taken from Wada et al. (2007) in which 
results of a deep 15 ^m survey during the performance 
verification phase are presented. In section 2, we de- 
scribe our data and reduction method using the IRC mon- 
itor observations. For detailed analysis of 18 /xm-selected 
galaxies, we focus on galaxies likely to lie at z ^ 0.5, i.e. 
the main targets of interests in the NEP survey. Then, 
we analyse these galaxies with the SED radiative trans- 
fer model of starbursts from Takagi et al. (2003b, a). In 
sections 3 and 4, we give our selection method for galax- 
ies at z ^ 0.5 and the SED fitting method, respectively. 
Results are described in section 5. In section 6, we at- 
tempt to identify silicate-break galaxies in our sample. 
We summarise our study in section 7. Throughout this 
paper, we adopt the cosmology of Qm — 0.3, f^A = 0.7 and 
Hq = 70kmsec~^Mpc^^. All magnitudes are given in the 
AB system, unless otherwise explicitly noted. 

2. Observations and data reduction 

2.1. Data from AKARI 

We have used the IRC monitor observations near 
the NEP with a field centre of a = 17h55'"24^ 
5 = -h66°37'32" (J2000), covering one field of view 
(~ 10' X 10') by all 9 IRC photometric bands, i.e. 
iV2, N3, N4, S7, S9W, 511, L15, L18W and L24. This 
field is regularly monitored with the IRC, in order to check 
the stability of the detectors. Wada et al. (2007) report 
a deep survey of this field at 3 and 7 /im, which was con- 
ducted during the performance verification phase. The ob- 
servations were done with the Astronomical Observation 
Template (ACT) three filter mode, IRC03 (see Onaka et 
al. 2007 for details). During one pointed observation, 
three IRC channels, NIR, MIR-S and MIR-L can operate 
simultaneously, covering different wavelength ranges. The 
NIR and MIR-S channels share the same field of view, 
while the MIR-L channel observes a region ~ 20' away 
from the field centre of NIR/MIR-S. In this paper, we use 
the data from 17 NIR/MIR-S pointed observations con- 
ducted between 19th June and 4th December 2006, and 5 
MIR-L pointed observations between 28th June and 30th 
July 2006. We summarise the observation parameters in 
Table 1. 

2.2. Reduction 

We reduced the data with the IRAF-based IRC pipeline 
software^ (Onaka et al. 2007), in order to remove and 
correct the basic instrumental effects: dark subtraction, 
linearity correction, distortion correction and flat field- 
ing. With the AOT IRC03, there is a sufficient number 
of frames per band (6-9 frames) during one pointed ob- 
servation to coadd the MIR-S data and also the MIR-L 
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data. We coadded the MIR-S frames from each pointed 
observation using the IRC pipeline in which sky positions 
of detected point sources arc matched to calculate the rel- 
ative astromctric shift between the frames. We did this 
after subtracting the sky background. For MIR-L, it was 
found that the sky background is too high and sources are 
too faint to match the positions of point sources. We used 
an optional command of the IRC pipeline (coaddLusingS) 
to coadd the MIR-L images, in which the relative astro- 
metric shifts of the MIR-L images are derived from those 
of the MIR-S images by adopting the calibration of the 
relative offset between the MIR-L and the MIR-S chan- 
nels. Thus we obtained coadded images of MIR-S and 
MIR-L for each pointed observation. 

In order to coadd images from different pointed obser- 
vations, we used the publicly available software SWarp'^, 
after putting the world coordinate system (WCS) on each 
image. We put the WCS for the MIR-S images using the 
IRC pipeline command putwcs in which 2MASS sources 
are identified and the astromctric solution is derived from 
them. For MIR-L, we manually identified three bright 
sources detected both in the MIR-L images and in the fi- 
nal Sll image, then determined the WCS for the MIR-L 
images. We used a median coaddition to produce the final 
images. 

For the NIR channel, we put the WCS into each indi- 
vidual frame using putwcs. This was successful for all of 
the NIR frames, although there is a large number of cos- 
mic ray hits in each frame. Again, SWarp was used to 
obtain the final images with median coaddition. 

2.3. Source detection and photometry on the IRC images 

Source detection was performed in a circular area of the 
final L18W image which has a radius of 4arcmin cover- 
ing the deepest region. We used SExtractor (Berlin & 
Arnouts, 1996) in which two or more connected > 3cr pix- 
els are identified as a source. This threshold results in 
the detection of sources with the signal-to-noise ratio of 
^ 4.5. This parameter setting for SExtractor is a result of 
our optimization to detect reliable point sources and avoid 
fainter extended sources which are likely to be residuals of 
the background subtraction. We have detected 72 sources 
in the area of 50.2 arcmin^. 

Since all of the 18 /im sources are found to have a coun- 
terpart at other IRC bands, the reliability of the detected 
sources is very high. We applied the same source detec- 
tion method for the negative L18W image, and detected 
only 5 sources. From an inspection of the negative image, 
it is obvious that they are not point sources like those 
detected in the positive image, but rather a result of the 
negative extended pattern caused by the residuals of the 
sky subtraction. We also used the negative L18W image 
to estimate the typical photometric error. On the neg- 
ative L18W image, we performed aperture photometry 
at random source-free positions (determined in the pos- 
itive image) using the same procedure as with detected 
sources (sec below). The distribution of flux measured 

^ http: / /terapix. iap. f r/rubrique .php?id_rubrique=49 
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at random positions has the standard deviation of 25 ^Jy, 
which represents a typical photometric error. We note 
that the faintest source in our catalogue has the flux of 
129±27/iJy. 

In the other IRC bands, we identified the counterparts 
of the detected 18 /Ltm sources by using the sky coordi- 
nates from SExtractor. First, we determined the centroid 
of the counterparts (with gcntrd.pro in the IDL astron- 
omy user's library^) using the sky coordinates of the 18 /im 
sources as the initial approximation of the centroid. If the 
newly-determined centroid is shifted from the original po- 
sition by > 3", we rejected such a source as the counterpart 
of the 18 ^m source. We found that this method results in 
a higher success rate of finding the counterparts of 18 /xm 
sources than the comparison of source catalogs from each 
IRC band with a certain search radius. We visually in- 
spected the result of the cross-identifications in postage 
stamp images (shown in Figures 3 &; 4 for the colour se- 
lected sample, see below), and verified the correctness of 
the cross-identifications. 

We performed aperture photometry at the centroid po- 
sitions of the 18 yum sources and their counterparts in our 
IRC images. We used aperture radii of 2 and 3 pixels 
for the NIR and the MIR-S/L images, respectively. Note 
that in the flux calibration of the IRC images with stan- 
dard stars, aperture radii of 10 and 7.5 pixels were used 
for the NIR and MIR-S/L channels, respectively. The 
aperture corrections used to estimate the total fluxes have 
been obtained from the empirical PSF derived from bright 
sources. For the NIR channel, we used bright sources in 
our images. For the MIR-S and the MIR-L channels, IRAS 
sources have been used which are detected in the obser- 
vations for the NEP-Wide survey. 

We found that the measured fluxes in the coaddcd NIR 
and MIR-S images are systematically lower than those be- 
fore coaddition. This is because we coadded the images, 
which have elongated PSFs, with various position angles 
using a simple median. Using the empirical PSFs, we gen- 
erated 'coadded' PSFs adopting the same position angle 
of the actual observations. The correction factors for this 
effect were estimated from aperture photometry of PSFs 
before and after coaddition. The correction factors are 
found to be 5% for N2, 7% for N3 and 7V4, and 3 % for 
the MIR-S channel. 

We have not applied colour corrections for the IRC pho- 
tometry at this initial stage of data reduction activity. 
The total error on the IRC photometry is estimated about 
10 - 20 %. 

2.4- Optical identification 

The IRC monitor fleld lies inside the NEP-Deep survey 
area and is covered by the deep BVRi'z' optical observa- 
tions with Subaru/S-cam, reaching B = 28.4 to z' = 26.2 
AB mag in 5a (Wada et al. 2007, in prep). This field 
is also covered by KPN02.1m/FLAMING0S at J- and 
A's-band with a depth of J = 21.6 and Kg = 19.9 (3tT) in 
Vega magnitudes (Ficld-NE in Imai et al., 2007). 

* http : //Idlastro .gsf c .nasa. gov/contents .html 



Since the astrometric accuracy is better for shorter 
wavelengths in the IRC bands, we adopted the posi- 
tions from the NIR channel to find the optical counter- 
parts of the 18 ^m sources. We used the catalogue of 
FLAMINGOS AT^-band sources with optical Subaru/S- 
cam identifications taken from Imai et al. (2007). We 
identified the nearest AT^-band source as the 18 /xm source 
counterpart. The resulting identifications were all visu- 
ally inspected. If the angular separation is larger than 3", 
wc reject the counterpart without visual inspection. We 
found no AT^-band counterparts for three 18 fim sources, 
for which we searched for the nearest optical counter- 
part from the Subaru/S-cam catalogue from Wada et al. 
(2007, in prep). A typical angular separation between 
the IRC positions and the optical positions is found to be 
0.6" ±0.2". 

3. Sample selection 

We have detected 72 sources in total at 18 /im, which 
include stars and galaxies at various redshifts. Since the 
main targets of this study and of the NEP surveys are star- 
forming galaxies at moderate - high redshifts, we adopt 
a simple colour cut to reduce the contamination from low 
redshift {z ^ 0.5) galaxies and stars. 

We did this selection by following Huang et al. (2004) 
in which z > 0.6 galaxies are selected using only the 
K — [3.6/im] colour by utilizing the 1.6 /im bump. Similar 
colour selection is possible with the IRC bands, i.e. A^2 — 
A^3 and/or A^2 — A^4. In order to confirm the effectiveness 
of this single colour selection, we derived photometric red- 
shifts of the 18 fim sources from the ground-based optical- 
NIR SEDs by using hyperz (Bolzonclla et al., 2000). We 
adopted the SED templates of Bruzual & Chariot (2003) 
distributed with the source code of hyperz. In Figure 1, 
we plot the derived photometric redshifts against A^2 — A^ 3 
and A^2 — NA colours. There is a clear correlation between 
Zhyperz and the AKARI NIR colours. This correlation con- 
firms that a rough selection of z ^ 0.5 galaxies is possible 
using only the single AKARI NIR colour. We prefer sin- 
gle colour cuts to photometric redshifts for the selection of 
galaxies at z ^ 0.5, due to their simplicity and availability 
over the entire field of the NEP survey. 

In the following we focus on the 25 galaxies satisfying 
the colour cut A^2 — A^3> 0.1 which are likely to lie at z > 
0.5. In this sample, we include ID65 which nevertheless 
has N2 — N3 < 0.1 with large errors, since its N2 — NA is 
similar to those of the other A^ 2 — A^3 selected galaxies and 
also its overall SED indicates that N2 flux is likely to be 
overestimated. As shown in Figure 1, a similar selection 
of galaxies is possible with A^2 — A^4 > —0.5. 

We further divided the 25 sources into two sub-classes, 
in order to separate AGN candidates from normal star- 
forming galaxies. With Spitzer surveys, the IRAC colour- 
colour plot of [3.6] — [5.8] vs [4.5] — [8.0] has been proposed 
for use in selecting AGN candidates (Lacy et al., 2004). 
Similarly, Oyabu et al. (2007) show that X-ray detected 
AGNs are found to be systematically red in both A^2 — A^ 4 
and A^3 — S7 indicating power-law SEDs. Based on their 
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result, we adopt a colour cut of iV3 — 57 > —0.2 for AGN 
candidates, which results in 11 AGN candidates in the 
18 urn sample with N2 - N3 > 0.1. In Figure 3 (4) we 
show postage stamp images of 14 (11) 7V2 — 7V3-selected 
star-forming galaxies (AGN candidates) at R- and Kg- 
band and 9 IRC bands. 

4. SED fitting method 

We analyse the SEDs of both 14 normal star-forming 
galaxies and 11 AGN candidates using the SED model 
of StarBUrsts with Radiative Transfer (SBURT; Takagi 
et al., 2003a, 2004). Since the contribution from AGN 
is not taken into account in SBURT, it is expected that 
the SEDs of AGN candidates show an excess in the MIR, 
owing to hot dust around the central massive black hole, 
while the best-fit model for normal star-forming galaxies 
would be largely acceptable. 

Wc find the best-fit SED model from the prepared set 
of 540 SED models at various redshifts with different star- 
burst age (0.01 - 0.6 Gyr), the compactness of the star- 
burst region 0.3 < 8 < 5.0 and the extinction curve (MW, 
LMC and SMC type)^. The compactness of starbursts is 
defined by r = 8(Mh./1O^M0) 2 [kpc], where r and Af* are 
the radius and stellar mass of the starburst region, respec- 
tively. Note that the attenuation due to dust or Ay is a 
function of starburst age, peaking at '^t/to where to is an 
evolutionary time-scale of starbursts taken to be 100 Myr. 
Also, the more compact the starburst region, i.e. smaller 
8, the larger the attenuation. The adopted SED models 
are the same as those used in Takagi et al. (2004) for the 
SED analysis of submm galaxies. 

Wc searched for the best-fit SED model by x^- 
minimization. Considering the flux uncertainty with the 
IRC images and systematic errors between ground-based 
and space photometry, we set minimum flux errors at each 
photometric band, i.e. 5% for optical bands {BVRi'z'), 
10% for J and K„ 15% for 7V2 and N3, 30% for NA and 
S7, and 20% for the other longer IRC bands. We have 
put larger minimum errors on iV4 and S7 bands, since a 
systematic error on calibration is indicated from the SED 
of stars in our field (see also Lee et al. 2007). Since these 
minimum errors are chosen as the most conservative limit, 
the derived values of can be regarded as the lower lim- 
its. 

At rest frame UV wavelengths, the SEDs of starbursts 
could depend on the small scale geometry of gas/dust and 
stars resulting in the leakage of UV photons and affect- 
ing the resulting SED (c.f. Takagi et al., 2003a). Also 
note that the extinction curve is not well understood in 
starbursts in the UV. Therefore, considering these uncer- 
tainties in the SED model, we quadratically added an ad- 
ditional 20% error for data at rest-frame UV wavelengths, 
<4000A, in order to reduce the statistical weight of data 
at the rest-frame UV. The adopted flux errors for the SED 
fitting arc shown in Figure 5 and 6. We reject the SED 

^ These SED models are available at 

http: //www. ir . Isas . jaxa. jp/'takagi/sedmodel/ . 
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model if the value of is significant at < 1% level. 
5. Results 

5.1. AGN /starburst diagnostics by the SED fitting 

Figures 5 and 6 show the best-fit SBURT model for 
normal star-forming galaxies and AGN candidates re- 
spectively, including those rejected with a large value 
(shown as dotted lines). As expected, most of the AGN 
candidates have power-law-like SEDs (see however ID25 as 
discussed below). According to the values, SBURT has 
no acceptable best-fit model for 4 out of 14 star-forming 
galaxies, while 7 out of 11 AGN candidates cannot be 
explained by SBURT. For the majority of galaxies with 
no acceptable SED model, the large is attributed to 
the excess of observed fluxes in more than 2 MIR bands, 
compared to the best-fit SED model. The detection of 
the MIR excess, probably due to the hot dust component 
around the central massive black hole, is reliable, thanks 
to the comprehensive wavelength coverage in the MIR. 
All the accepted models for AGN candidates (ID14, 35, 
64, and 72) underestimate the MIR fiuxes as well, but this 
does not much contribute to the x^ value, due to the large 
minimum error (30%) for A^4 and 5*7. Thus, if more ac- 
curate photometry becomes possible for these bands, the 
rejection rate for AGN candidates could reach ~100%. 
This result indicates that the simple colour selection of 
AGN candidates with AKARI bands is consistent with 
the more sophisticated diagnostics using the SED fitting. 

Assuming all of the AGN candidates selected with 
N2-m> 0.1 and Ni-S7 > -0.2 are indeed AGN, 
we estimate the AGN fraction of 18/im-selected galax- 
ies to be 15±4% (11/72) in number or ~12% in total 
fiux. Although the completeness correction is not yet ap- 
plied, this fraction is consistent with previous MIR surveys 
(e.g. Fadda et al., 2002; Caputi et al., 2006). For exam- 
ple, Fadda et al. (2002) estimate that AGNs contribute 
15±5 % of the total 15 /J,m flux in the Lockman Hole with 
a survey depth of 0.3 mJy at 15 /xm. Manners et al. (2004) 
similarly find an AGN fraction of ~ 19% from Chandra fol- 
low up of 0.8 - 6mJy 15 /im sources from the European 
Large Area ISO survey. 

5.2. Photometric redshift using MIR features 

MIR photometry, such as Spitzer 24 /im band data, are 
rarely used to derive photometric redshifts (see however 
Perez-Gonzalez et al., 2005). By using the superb MIR 
coverage with the IRC photometry, it may be possible to 
constrain the redshifts of distant infrared galaxies by using 
the PAH and silicate absorption features. In Figure 5, we 
find that the steep rise in fiux corresponding to the blue 
side of the 6.2 /im feature is successfully traced by the 
well-sampled IRC photometric bands; the best example 
of this can be seen in IDIO. This is not possible with 
Spitzer photometry, which has a large gap between 8 and 
24 //m. This new capability of the IRC may be useful for 
determining photometric redshifts for optically faint and 
featureless galaxies, such as young and dusty starburst 
galaxies with less prominent 4000 A breaks. However, it is 
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still necessary to calibrate the MIR photometric redshifts 
using spectroscopic follow-up observations. 

In Figure 7, we show the comparison of photomet- 
ric redshifts derived with hyperz using the ground-based 
optical-NIR photometry and those with the optical-to- 
IRC bands using SBURT. The majority of these photo- 
metric redshifts are consistent with one another within 
the 99% confidence limits. Significant discrepancies can 
be seen in ID14, 24, and 72. Among these, ID14 and 72 
are AGN candidates, and ID24 is a star- forming galaxy. 
For AGN candidates, photometric redshifts from both hy- 
perz and SBURT would not be reliable, since the adopted 
models do not include an AGN component. The star- 
forming galaxy ID24 has a nearly power-law-like SED in 
the optical-NIR bands, and therefore the photometric red- 
shift using only the optical-NIR SED is likely to be un- 
reliable. Photometric redshifts with the IRC bands could 
be most useful for this type of galaxy. 

As suggested by Takagi & Pearson (2005), the PAH and 
silicate absorption features of galaxies at z 1 could be 
identified using IRC 15, 18 and 24 /^m photometry. Due 
to the relatively shallow depth of our 24 /.utl image, their 
'silicate-break' selection technique is only partly applica- 
ble to our data. We discuss the application of this tech- 
nique for our sample in section 6. 

5.3. The nature of 18 fim-selected galaxies 

Here we briefly describe the physical properties of the 
18 /xm-selected galaxies with N2 — N3 > 0.1 based on the 
results of the SED fitting using the SBURT model. More 
detailed analyses of AKARI MIR-selected galaxies will be 
given elsewhere. 

The average 18/im flux of star-forming galaxies with 
N3- S7 < -0.2 is 240 ^Jy. In Table 2, we summarise 
the derived quantities from the SED fitting of the galax- 
ies with acceptable (i.e. not significant at < 1% level), 
along with the best-fit parameters. The average redshift 
of these galaxies is (z) = 1.0. Infrared luminosities range 
from 10^^ to 3 X lO^^L©, and therefore they are classi- 
fied as LIRGs or ULIRGs. From the SBURT model, a 
typical stellar mass of these (U)LIRGs is derived to be 
3 X lO^°Af0, including galaxies as massive as ^ IO^^Mq. 

It is interesting to note that no best-fit SED model 
shows prominent silicate absorption feature, i.e. they 
are less obscured objects, unlike nearby ULIRGs. There 
may be a systematic difference in the dust attenuation 
between low- and high-z objects with the same luminos- 
ity. At high redshifts, strong silicate absorption similar 
to nearby ULIRGs (e.g. Armus et al., 2007) is found 
in more luminous objects such as submm galaxies or 
MIR-selected galaxies with > IO^^Lq (e.g. Houck et al., 
2005). A systematic difference in the SEDs between 
ULIRGs at z < 1 and submm galaxies (i.e. Hy/ULIRGs 
at z ^ 2) is already reported as a systematic difference 
in dust temperatures (Chapman et al., 2005; Yang et al., 
2007) and in the UV-submm SEDs (Takagi et al., 2004). 
Also, Reddy et al. (2006) claim that at a given bolometric 
luminosity the obscuration of galaxies is generally ^ 10 
times smaller at z ~ 2 than at z ~ 0. The other possibility 



is that our results indicate a selection effect, since 
ULIRGs with strong PAH emission should be brighter 
than obscured ULIRGs. A systematic difference in the 
SED of galaxies at different redshifts could be one of the 
most important topics for further work in the NEP survey. 

Cases worthy of comment 

IDIO: The PAH emission of this galaxy is well sampled 
with the IRC photometry given the photometric redshift 
of zphot = O-Gto.l- With MIR-L bands, the PAH 11.2 /^m 
feature is covered. From the SED fitting, wc find that the 
resulting is 24.9 and half of this value is explained 
by the L15 and L18W bands alone, since the model over- 
estimates the observed fluxes. This may indicate that the 
PAH emission of this galaxy has anomalous inter-band 
strength ratios. Certainly, we can put more strong con- 
straints on this inter-band strength ratio once the spectro- 
scopic redshift becomes available. With the mapping of 
the nearby spiral galaxy NGC6946 using AKARI/IRC at 
7 and 11 jim and also the MIR spectroscopy, Sakon et al. 
(2007) report a strong variation of PAH 7.7/11.2 /im ra- 
tio, depending on the physical properties of interstellar 
medium. In the distant universe, galaxies with anoma- 
lous PAH inter-band ratio may exist in greater number 
than expected from the local universe. 

ID25: Although this galaxy is classified as an AGN 
candidate, its SED is not characterized by a simple power- 
law. The 1.6/im bump falls in between A^2, TVS and A^4 
bands, suggesting that stellar emission dominates in this 
wavelength range. From 7/im, there is a very steep rise in 
flux towards longer wavelengths which cannot be repro- 
duced with the SBURT models. 

ID26: This is an extremely red object with i? — A' = 
4.6 (6.3 in Vega), and too red for the SBURT models 
within the reasonable parameter range for the SED fltting. 
Interestingly, the SED fltting with hyperz indicates that 
this is a young (40 Myr) and heavily obscured {Ay = 5) 
galaxy at z = 1. Also, this galaxy has the bluest L15 — 
L18W colour in our sample, i.e. good candidate of silicate- 
break galaxy at z ~ 1 (see section 6). 

ID65: This galaxy has a very steep power-law SED 
with the optical-to-IR slope of a ~ —2.5 when oc 
This is very close to the maximum limit on the steepness 
for AGN SEDs (e.g. Alonso-Herrero et al., 2006). 

6. Discussion 

Takagi & Pearson (2005) proposed a new selection tech- 
nique of high-z ULIRGs by using the AKARI/IRC multi- 
wavelength survey, in which the colour anomaly in the 
MIR bands due to the silicate absorption feature is used. 
Galaxies found with this method are called silicate-break 
galaxies. Kasliwal et al. (2005) adopted this technique 
by using Spitzer 16 /^m (IRS) and 24 /im (MIPS) obser- 
vations. They found 36 candidate silicate-break galaxies 
probably at 1 ^ z < 1.8 in 0.0392 deg^, corresponding ^ 920 
sources deg^. 

Here wc attempt to identify silicate-break galaxies with 
our data. As pointed out by Takagi & Pearson (2005), 



No. ] 

the key wavelength for this selection method is 24 /im. 
Unfortunately, our 24/xin image is not deep enough to em- 
ploy the silicate-break technique. Instead, we use Lib and 
L18W. 

In Figure 8, we show the colour-colour plot of N2 — N3 
vs. L15-L18W. According to Takagi & Pearson (2005), 
we could select ULIRGs at z 1 with L15 - L18W < 0, 
but we would suffer from the contamination from normal 
spiral galaxies at low redshifts. We use a colour cut of 
N2 — N3 > 0.1 again to exclude the contamination from 
low-z galaxies. 

We find that seven 18 /xm-sclccted galaxies satisfy N2 — 
iV3 > 0.1 and Lib - L18W < 0. Five out of seven, ID15, 
16, 26, 58 and 60, have iV3 - 5*7 < -0.2, and arc hkcly to 
be star-forming galaxies. Except for ID16 and 58, photo- 
metric redshifts (from both SBURT and hyperz) are con- 
sistent with the rcdshift expected from the silicate-break 
selection, i.e. z ~ 1. ID16 and 58 have Zphot = 0.6 and 
0.7, respectively. These low-z contaminants are expected, 
since the single colour cut N2 — N3 > 0.1 is only useful for 
a rough selection of galaxies at z ^ 0.5. Such contamina- 
tions would be reduced if we use a colour selection with 
18 and 24 //m as suggested by Takagi & Pearson (2005). 

As noted in section 5.3, no best-fit SED models for 
silicate-break galaxies have strong silicate absorption. So, 
they arc not silicate-break but 'PAH'-brcak galaxies. On 
the other hand, Takagi & Pearson (2005) show that the 
SED of a silicate-break galaxy selected with the Spitzer 
IRS peak-up imager can be reproduced by an SED model 
with a strong silicate absorption feature. This galaxy is 
a submm galaxy and therefore more luminous than our 
18 /zm-selected galaxies. With deeper and wider MIR-L 
observations, we would find ULIRGs at z ~ 1 with silicate 
absorption, i.e. true silicatc-brcak galaxies. More stud- 
ies on the nature of silicate-break galaxies will be given 
elsewhere using the main NEP survey. 

7. Summary 

We investigated AKARI 18 /.tm-selected sources by us- 
ing all nine IRC photometric bands, as a pilot study for 
the main AKARI NEP survey. By using the IRC observa- 
tions for the monitor field, we detected 72 18 /im sources 
from the area of 50.2 arcmin^. In order to select galax- 
ies at z ^ 0.5, i.e. the main targets for the NEP survey, 
we employed a simple colour selection, A^2 — A^3 > 0.1 
and obtained a sub-sample of 25. AGN candidates in 
this sub-sample are identified with another colour cut, 
N3~S7 > —0.2, which implies a hot dust component. The 
resulting sample has 14 star-forming galaxies at z ^ 0.5 
and 11 AGN candidates with power-law like SEDs. 

We analysed both star-forming galaxies and AGN can- 
didates with a radiative transfer SED model of starbursts, 
SBURT (Takagi et al., 2003a). Since the SBURT model 
does not include an AGN component, the presence of an 
AGN can be assessed by the excess of observed flux at the 
MIR wavelengths due to hot dust around the central mas- 
sive black hole. We found that our AGN candidates show 
MIR excesses compared to the best-fit SED models, most 
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of which are rejected at the significance of < 1% by the re- 
sulting large value. This means that the simple colour 
selection for AGN candidates is consistent with the more 
sophisticated diagnostics with the SED model. An inter- 
esting future project would be SED analysis with more 
complicated models having both starbursts and AGN. 

With the SED fitting, we obtained photometric red- 
shifts using optical-to-MIR SEDs, which are largely con- 
sistent with ground-based optical-NIR photometric red- 
shifts. We demonstrated that the steep rises in flux cor- 
responding to the blue-side of the PAH 6.2 /tm feature 
and also the flux dip between the PAH 7.7 and 11.2 /tm 
features are recognizable in the resulting SEDs by IRC. 
The IRC all-band photometry is probed to be useful to 
constrain the redshift of infrared galaxies using the MIR 
spectral features. 

Typically, A^2 — A^3-sclected 18 /im sources are 
(U)LIRGs at z - 1 with masses ~ 3 x lO^^M©. No 
ULIRG in our sample has strong silicate absorption 
according to the best-fit SED model. This may suggest 
a systematic difference in the nature of ULIRGs at low 
and high redshifts, or may be a simple selection effect, 
owing to the faintness of ULIRGs with strong silicate 
absorption. We could test these hypotheses with deeper 
IRC all-band surveys such as the NEP-Deep survey. 

This pilot study demonstrates well the interesting new 
capabilities of the AKARI NEP survey, compared to the 
existing infrared surveys by Spitzer. The IRC all-band 
photometry would be a powerful and unique tool to con- 
strain the redshifts of optically faint heavily obscured 
galaxies. We may find many candidate examples of in- 
frared galaxies with anomalous PAH inter-band strength 
ratio. Furthermore, unique samples of ULIRGs at z ~ 1 
will be constructed with the silicate-break selection tech- 
nique. Such samples could play an important role in re- 
vealing the nature of high-z ULIRGs, which should pro- 
vide vital clues on the process of galaxy formation. 
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Table 1. Summary of observations 
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Table 2. Results of the SED fitting with SBURT") 
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Fig. 1. Photometric redshifts derived from the ground-based optical-NIR SED by using hyperz as a function of the IRC NIR colour, 
N2 - N3 and N2 - Ni. Errors of 

^hyperz the 99% Confidence limit. Encircled dots are for galaxies with — N3 > 0.1. Most of 
the galaxies with — N3 > 0.1 can be selected with — > —0.5. Large crosses indicate that the best-fit SED model is rejected 
with the significance of < 1%. The tight correlation between photometric redshifts and NIR colours shows that a NIR colour cut 
can closely mimic a cut with the photometric redshift. 




Fig. 2. N2 — N3 vs. — S7 colour-colour diagram to classify 18 /im-sclcctcd galaxies into low-z (z < 0.5) galaxies, high-z star- 
forming galaxies and AGN candidates. 
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Fig. 3. Postage stamp image of 18 /xm-selected galaxies with N2 — N3 > 0.1 and A'S — 57 < —0.2, i.e. star-forming galaxies at 
z > 0.5. From left to right, we show R, Ks, N2, N3, NA, S7,S9W, 511, L15, L18W and L24 band images. Identified L18W sources 
are indicated with black and white circles with the radii of 3.5" and 3" respectively. The north is up and the east is to the left. 




Fig. 4. Same as Figure 3, but for galaxies with N2~ N3> 0.1 and A^S - 57 > -0.2, i.e. AGN candidates. 
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Fig. 5. Results of SED fitting for 18 ^im-selected galaxies with N2 - N3 > 0.1 and A^S - 57 < - 
z ^ 0.5. The best-fit SED models rejected with the significance of < 1% are shown in dotted lines. 
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Fig. 6. Same as Figure 5, but for galaxies with A^2 - TVS > 0.1 and A^3 - 57 > -0.2, i.e. AGN candidates. 
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Fig. 7. Comparison of photometric redshifts derived from hyperz and SBURT. Errors on photometric redshifts are in 99% confidence 
Hmit. The numbers indicate the object ID of corresponding data points. Squares indicate star-forming galaxies. 
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Fig. 8. Colour-colour plot to identify silicate-break galaxies. Squares indicate star-forming galaxies which satisfy selection criteria 
of silicate-break galaxies. The numbers indicate the object ID of corresponding data points. 



